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Abstract. The phenomenon of crossed Andreev reflection (CAR) is known to play a 
key role in non-local electron transport across three-terminal normal-superconducting- 
normal (NSN) devices. Here we review our general theory of non-local charge trans- 
port in three-terminal disordered ferromagnet-superconductor-ferromagnet (FSF) 
structures. We demonstrate that CAR is highly sensitive to electron spins and 
yields a rich variety of properties of non-local conductance which we describe non- 
perturbatively at arbitrary voltages, temperature, degree of disorder, spin-dependent 
interface transmissions and their polarizations. We demonstrate that magnetic ef- 
fects have different implications: While strong exchange field suppresses disorder- 
induced electron interference in ferromagnetic electrodes, spin-sensitive electron 
scattering at SF interfaces can drive the total non-local conductance negative at 
sufficiently low energies. At higher energies magnetic effects become less important 
and the non-local resistance behaves similarly to the non-magnetic case. Our results 
can be applied to multi-terminal hybrid structures with normal, ferromagnetic and 
half-metallic electrodes and can be directly tested in future experiments. 



1.1 Introduction 

In hybrid NS structures quasiparticle current flovifing in a normal metal is 
converted into that of Cooper pairs inside a superconductor. For quasiparticle 
energies above the superconducting gap e > A this conversion is accompanied 
by electron- hole imbalance which relaxes deep inside a superconductor. At 
subgap energies e < A the physical picture is entirely different. In this case 
quasiparticIe-to-Cooper-pair current conversion is provided by the mechanism 
of Andreev reflection (AR) [1]: A quasiparticle enters the superconductor from 
the normal metal at a length of order of the superconducting coherence length 
,fs, forms a Cooper pair together with another quasiparticle, while a hole goes 
back into the normal metal. As a result, the net charge 2e is transferred 
through the NS interface which acquires non-zero subgap conductance down 
to r = o [2]. 
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AR remains essentially a local effect provided there exists only one NS 
interface in the system or, else, if the distance between different NS inter- 
faces greatly exceeds the superconducting coherence length ^. If, however, the 
distance L between two adjacent NS interfaces (i.e. the superconductor size) 
is smaller than (or comparable with) ^, two additional non-local processes 
come into play (see Fig. 1.1). One such process corresponds to direct electron 
transfer between two N-metals through a superconductor. Another process is 
the so-called crossed Andreev reflection [3, 4] (CAR): An electron penetrating 
into the superconductor from the first N-terminal may form a Cooper pair 
together with another electron from the second N-terminal. In this case a 
hole will go into the second N-mctal and AR becomes a non-local effect. This 
phenomenon of CAR enables direct experimental demonstration of entangle- 
ment between electrons in spatially separated N-electrodes and can strongly 
influence non-local transport of electrons in hybrid NSN systems. 




Fig. 1.1. Two elementary processes contributing to non-local conductance of an 
NSN device: (1) direct electron transfer and (2) crossed Andreev reflection. 



Non-local electron transport in the presence of CAR was recently in- 
vestigated both experimentally [5, 6, 7, 8, 9, 10, 11, 12] and theoretically 
[13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26] (see also further ref- 
erences therein) demonstrating a rich variety of physical processes involved 
in the problem. It was shown [13] that in the lowest order in the interface 
transmission and at T = CAR contribution to cross-terminal conductance 
is exactly canceled by that from elastic electron cotunneling (EC), i.e. the 
non-local conductance vanishes in this limit. Taking into account higher order 
processes in barrier transmissions eliminates this feature and yields non-zero 
values of cross-conductance [14]. 

Another interesting issue is the effect of disorder. It is well known that 
disorder enhances interference effects and, hence, can strongly modify local 
subgap conductance of NS interfaces in the low energy limit [27, 28, 29]. Non- 
local conductance of multi-terminal hybrid NSN structures in the presence of 
disorder was studied in a number of papers [15, 16, 17, 18, 24, 25]. A general 
quasiclassical theory was constructed by Golubev and the present authors [24] . 
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It was demonstrated that an interplay between CAR, quantum interference of 
electrons and non-local charge imbalance dominates the behavior of diffiisive 
NSN systems being essential for quantitative interpretation of a number of 
experimental observations [7, 9, 10]. In particular, strong enhancement of 
non-local spectral conductance was predicted at low energies due to quantiim 
interference of electrons in disordered N-terminals. At the same time, non- 
local resistance R12 remains smooth at small energies and, furthermore, was 
found to depend neither on parameters of NS interfaces nor on those of N- 
terminals. At higher temperatures R12 was shown to exhibit a peak caused 
by the trade-off between charge imbalance and Andreev reflection. 

Yet another interesting subject is an interplay between CAR and Coulomb 
interaction. The eflfect of electron-electron interactions on AR was investigated 
in a number of papers [29, 30, 31]. Interactions should also affect CAR, e.g., by 
lifting the exact cancellation of EC and CAR contributions [19] already in the 
lowest order in tunneling. A similar effect can occur in the presence of external 
ac flelds [23]. A general theory of non-local transport in NSN systems with 
disorder and electron-electron interactions was very recently developed by 
Golubev and one of the present authors [26] direct relation between Coulomb 
effects and non-local shot noise. In the tunneling limit non-local differential 
conductance is found to have an S-like shape and can turn negative at non- 
zero bias. At high transmissions CAR turned out to be responsible both for 
positive noise cross-correlations and for Coulomb anti-blockade of non-local 
electron transport. 

An important property of both AR and CAR is that these processes should 
be sensitive to magnetic properties of normal electrodes because these pro- 
cesses essentially depend on spins of scattered electrons. One possible way to 
demonstrate spin-resolved CAR is to use ferromagnets (F) instead of normal 
electrodes ferromagnet-superconductor-ferromagnet (FSF) structures. First 
experiments on such FSF structures [5] illustrated this point by demonstrat- 
ing the dependence of non-local conductance on the polarization of ferro- 
magnetic terminals. Hence, for better understanding of non-local effects in 
multi-terminal hybrid proximity structures it is necessary to construct a the- 
ory of spin-resolved non-local transport. In the case of ballistic systems in the 
lowest order order in tunneling this task was accomplished in [13]. Here we 
will generalize our quasiclassical approach [14, 24] to explicitly focus on spin 
effects and construct a theory of non-local electron transport in both ballis- 
tic and diffusive NSN and FSF structures with spin-active interfaces beyond 
lowest order perturbation theory in their transmissions. 

The structure of the paper is as follows. In Sec. 2 we will describe non-local 
spin-resolved electron transport in ballistic NSN structures with spin-active 
interfaces. In Sec. 3 we will further extend our formalism and evaluate both 
local and non-local conductances in SFS structures in the presence of disorder. 
Our main conclusions are briefly summarized in Sec. 4. 
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1.2 Spin-resolved transport in ballistic systems 

Let us consider three-terminal NSN structure depicted in Fig. 1.2. We will 
assume that all three metallic electrodes are non-magnetic and ballistic, i.e. 
the electron clastic mean free path in each metal is larger than any other 
relevant size scale. In order to resolve spin-dependent effects we will assume 
that both NS interfaces are spin-active, i.e. we will distinguish "spin-up" and 
"spin-down" transmissions of the first (Di-^ and Di^) and the second (-02-^ and 
D21) SN interface. All these four transmissions may take any value from zero 
to one. We also introduce the angle between polarizations of two interfaces 
which can take any value between and 27r. 

In what follows effective cross-sections of the two interfaces will be denoted 
respectively as Ai and A2- The distance between these interfaces L as well as 
other geometric parameters arc assumed to be much larger than \J~A\j.. i.e. 
effectively both contacts are metallic constrictions. In this case the voltage 
drops only across SN interfaces and not inside large metallic electrodes. 



For convenience, we will set the electric potential of the S-electrode equal to 

zero, = 0. In the presence of bias voltages Vi and V2 applied to two normal 
electrodes (see Fig. 1.2) the currents I\ and I-i will flow through SNi and SN2 
interfaces. These currents can be evaluated with the aid of the quasiclassical 
formalism of nonequilibrium Green-Eilenberger-Keldysh functions [32] g^'^'^ 
which we briefly specify below. 

1.2.1 Quasiclassical equations 

In the ballistic limit the corresponding Eilenberger equations take the form 




Fig. 1.2. Schematics of our NSN device. 



eT3 + eV{r,t)- A{r,t),g 



+ivF'Vg 



'{PF,e,r,t) = 0, 



(1.1) 
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where [a, b] = ab — ba, e is the quasiparticle energy, pp = mvp is the electron 
Fermi momentum vector and is the Pauh matrix in Nambu space. The 
functions g^'^'^ also obey the normalization conditions (g^)^ = {g^)^ = 1 
and g^g^ + g^g^ = 0. Here and below the product of matrices is defined as 
time convolution. 

Green functions g^'^'^ and A are 4x4 matrices in Nambu and spin spaces. 
In Nambu space they can be parameterized as 

/„R,A,K fR,A,K\ / n 



9 ' ' - yR,A,K gR,A,K J . ^ - ) ' ^^-^^ 

where g^'^'^ , jR,a,k ^ jR,A,K ^ gR,A,K g^j^g 2x2 matrices in the spin space, A 
is the BCS order parameter and CTj are Pauli matrices. For simplicity we will 
only consider the case of spin-singlet isotropic pairing in the superconducting 
electrode. The current density is related to the Keldysh function g^ according 
to the standard relation 



3{'r,t) = ^^ j de{vFSp[f3g'^{pF,e,r,t)]) , 



(1.3) 



where; .¥0 = mpp/2n'^ is the density of state at the Fermi level and angular 
brackets (...) denote averaging over the Fermi momentum. 



1.2.2 Riccati parameterization 

The above matrix Green- Keldysh functions can be conveniently parameterized 
by four Riccati amplitudes 7^'"*, ^^-^ and two "distribution functions" , 
(here and below we chose to follow the notations [33]): 

9 ='^N ( ^k'_^r^k^a ]N , (1.4) 



where functions j^'^ and 7^'^ are Riccati amplitudes 

(1 _i_ ^/-R,'4^-,-R,A n R,A 

'%«7a _,'^^R,A,nA] (1-5) 

and N^'^ are the following matrices 



^ -[ ^i_;yR,AjR,Ayl)- 

With the aid of the above parameterization one can identically transform the 

quasiclassical equations (1.1) into the following set of effectively decoupled 
equations for Riccati amplitudes and distribution functions [33] 
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ivpV-i^^^ + [e + eV{r,t)]^^'^ + ^^^^[e - eV{r,t)] 

= 7«'^Zi*ia27-^'^ - Aia2, 
ivpS/j^'^ - [£ - ey(r,i)]7^'^ - ^^'^[e + eV{r,t)] 



(1.7) 



(1.8) 



ivpVx^ + [e + eV{r, t)]x^ 



.K 



x^[e + eV{r,t)] 
x'^Aia2^^ = 0, 



,K 



(1.9) 



ivpVx'^ - [£ - eV{r, t)]i^ + [e - eV{r, t)] 



(1.10) 



Depending on the particular trajectory it is also convenient to introduce a 
"replica" of both Riccati amplitudes and distribution functions which - again 
following the notations [33, 34] - will be denoted by capital letters F and 
X. These "capital" Riccati amplitudes and distribution functions obey the 
same equations (1.7)-(1.10) with the replacement 7 — ?■ and x ^ X. The 
distinction between different Riccati amplitudes and distribution functions 
will be made explicit below. 

1.2.3 Boundary conditions 




Fig. 1.3. Riccati amplitudes for incoming and outgoing trajectories from the both 
sides of the interface. 

Quasiclassical equations should be supplemented by appropriate bound- 
ary conditions at metallic interfaces. In the case of specularly reflecting spin- 
degenerate interfaces these conditions were derived by Zaitsev [35] and later 
generalized to spin-active interfaces [36] , see also [37] for recent review on this 
subject. 

Before specifying these conditions it is important to emphasize that the ap- 
plicability of the Eilenberger quasiclassical formalism with appropriate bound- 
ary conditions to hybrid structures with two or more barriers is, in general, a 
non- trivial issue [38, 39]. Electrons scattered at different barriers interfere and 
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form bound states (resonances) which cannot be correctly described within 
the quasiclassical formalism employing Zaitscv boimdary conditions or their 
direct generalization. Here we avoid this problem by choosing the appropriate 
geometry of our NSN device, see Fig. 1.2. In our system any relevant trajec- 
tory reaches each NS interface only once whereas the probability of multiple 
reflections at both interfaces is small in the parameter AiA2/L'^ "C 1. Hence, 
resonances formed by multiply reflected electron waves can be neglected, and 
our formalism remains adequate for the problem in question. 

It will be convenient for us to formulate the boundary conditions directly 
in terms of Riccati amplitudes and the distribution functions. Let us consider 
the first NS interface and explicitly specify the relations between Riccati am- 
plitudes and distribution functions for incoming and outgoing trajectories, see 
Fig. 1.3. The boundary conditions for r{^, and Xi can be written in the 
form [34] 

r^ = rh^St,+tZl^,S+,„ (1.11) 

= S.^lMr + S,y^vit, (1-12) 
= r^^xfft + tM,it - a^ixfrat. (1.13) 

Here we deflned the transmission (t), reflection (r), and branch-conversion (a) 
amplitudes as: 

r^i = +[{Pvi)-'S+ - {l3?vr'S+,]-\l3?,r\ (1.14) 

til = -[Wvi)-'S+ - W^r)-'S+,]-'iPvi')-\ (1-15) 

ft = +(/3ri)-M^ii(/3ri)-' - S,rif3fn')-']-\ (1-16) 

it = -Wi'i')-'[SiiW^^)-' - S,Af3^'i')-T\ (1-17) 

a^^ = (r«5ii-5n7f)(;3«,)-\ (1-18) 

at = (^n')-'(^liA^ - 7i^^i+i), (1-19) 

where 

^5 = 4 - 7f ^+-7f , = Sj, - 7f 5,,7f , (1.20) 

Pfj = Sij - itS^^lf, Ptj = S% - itself. (1.21) 

Similarly, the boundary conditions for , F^ , and take the form: 

r^ = f^ii^s^^ + i^a^Sv,, (1.22) 

F^ = SWrt + St,,lt'tL (1-23) 
If = f^^if rf, + t^ii'^^tt - ~a'^ix^>at, (1-24) 



where 
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(1.25) 


til — 




(1.26) 


- 

' Ir — 




(1.27) 


- 

•'Ir — 


-{M'v)-\stSv,)-' - st,Svv)-']-\ 


(1.28) 




af, = (A^5+-5+7f)(/3fi,)-\ 


(1.29) 






(1.30) 



Boundary conditions for Fy'"^, P^'^, Xy and can be obtained from the 
above equations simply by replacing 1 O 1'. 

The matrices Sn, Sw , Svi, and Sw constitute the components of the 
<S-matrix describing electron scattering at the first interface: 




Fig. 1.4. Riccati amplitudes for incoming and outgoing trajectories for an NSN 
structure with two barriers. The arrows define quasiparticle momentum directions. 

We also indicate relevant Riccati amplitudes and distribution functions parameter- 
izing the Green-Keldysh function for the corresponding trajectory. 



In our three terminal geometry nonlocal conductance arises only from 
trajectories that cross both interfaces, as illustrated in Fig. 1.4. Accordingly, 
the above boundary conditions should be employed at both NS interfaces. 

Finally, one needs to specify the asymptotic boundary conditions far from 
NS interfaces. Deep in metallic electrodes we have 
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7f = 7f = it = it = 0, (1-32) 

= ho{e + eVx), x^ = -ho{£ - eV{), (1.33) 

71^ = 72'' = 72^ =72^ = 0, (1.34) 

a;f = /io(£ + eV2), = -/io(e - 6^2), (1-35) 

where /io(£) = tanh(£/2T) - equilibrium distribution function. In the bulk of 
superconducting electrode we have 

7^ = -a(£)i(T2, 7j1 = a*(£)i(T2, (1.36) 

xf,=-[l-\a{e)\'']ho{e), (1.37) 

7^5 = a(£)io-2, 7^ = -a*(£)i(72, (1.38) 

a;|^ = [1 - |a(£)fl/io(£), (1.39) 

where we denoted a(£) = — (e — Ve^ — /4^)/^. 
1.2.4 Green functions 

With the aid of the above equations and boundary conditions it is straight- 
forward to evaluate the quasiclassical Green-Keldysh functions for our three- 
terminal device along any trajectory of interest. For instance, from the bound- 
ary conditions at the second interface we find 

= w(£)^2fT2, (1.40) 

where A2 = S'2'2'0'2j22'2'''^2- Integrating Eq. (1.7) along the trajectory connect- 
ing both interfaces and using Eq. (1.40) as the initial condition we immediately 
evaluate the Riccati amplitude at the first interface: 

- ' l-iaA2A + e)Q ^^'^^^ 
Q^^_^f^^^mM^ n=^7^^. (1.42) 

Employing the boundary conditions again we obtain 

= iSivK-^ M2 + (0^26 + A)Q\ (725+ , , (1.43) 

f« = -iaSX,x02SvvK2t [1 - {aA2A + e)Q\ S-\,Svx, (1.44) 

where 

= (1 - a^A^Aj) - [£(1 + a^A,Aj) + Aa{Ai + A^)\ Q, (1.45) 

Ai=a2S+^,a2Svv. (1.46) 
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We also note that the relation (r^'^)+ = r^'^ and (7^'^)+ = 7^'^ makes 
it unnecessary (while redundant) to separately calculate the advanced Riccati 
amplitudes. 

Let us now evaluate the distribution functions at both interfaces. With 
the aid of the boundary conditions at the second interface we obtain 

X§ = S2'2'S+2, (1 - |a|2) ho{s) + S^'iS+^x^ 

— |ci| S'2'2'f2i5^/ ^22' a2S+,2,x^. (1.47) 

Integrating Eq. (1.9) along the trajectory connecting both interfaces with 
initial condition for , we arrive at the expression for Xy 

x^, =[l-{aA2A + e)Q]-^ X§ 

X (1 - tanh^ iL[2/vF) [1 - {aAiA + £)Q]+"' . (1.48) 

Then we can find distribution functions at the first interface. On the normal 
metal side of the interface we find 

= ril^f + t^ixMi + afzaff (l - \a\^) K{s) (1.49) 



where 



r^i = Su'K^i (1 - iaA2A + e)Q)S+,,S+; 



(1.50) 

-0(0^2 + {aA2e + A)Q)a2S+,ya2S+-^ , 

t^i = SivK-^\l - iaA2A + e)Q), (1.51) 
a^i = iSivK2j'{aA2 + (a^2£ + A)Q)a2S+^,. (1.52) 
The corresponding expression for is obtained analogously. We get 

= r^^^f fr - i^ii^i" (1 - ho{e) - a«xf^4". (1.53) 

where 



(1.54) 



S-,\S^,^,a2il-{aA2A + e)Q) 

-Sy\a2Syva{aA2 + (a^2e + A)Q) K^2a2Sty, 

ill = St,iS+,l!(T2{l - {aA2A + e)Q)K^^(T2S+^,, (1.55) 

ail = iaS+^(j2SvvK^^{l - (0^2^ + £)Q)- (1-56) 

Combining the above results for the Riccati amplitudes and the distribution 
functions we can easily evaluate the Keldysh Green function at the first in- 
terface. For instance, for the trajectory 1°"* (see Fig. 1.4) we obtain 

5i^„*=2(Xf -r«xfr«^), =2if. (1.57) 

The Keldysh Green function for the trajectory 1™ is evaluated analogously, 
and we get 

= 2xf , = 2(Xf - r«^f n^")- (1-58) 
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1.2.5 Nonlocal conductance: General results 

Now we are ready to evaluate the current /i across the first interface. This 
current takes the form: 



h = ir " (V^i) - ^ / d£ Sp(f35fo.* - f35f»), (1-59) 



where 



Go = ^^:2r2 (1-60) 



87i72M.A/'2 

is the normal state conductance of our device at fully transparent inter- 
faces, Pf7i(2) is normal to the first (second) interface component of the 
Fermi momentum for electrons propagating straight between the interfaces, 
A/i.2 = PpAi.2/4:TT define the number of conducting channels of the corre- 
sponding interface, Rg = 27r/e^ is the quantum resistance unit. 

Here /^^^(Fi) stands for the contribution to the current through the first 
interface coming from trajectories that never cross the second interface. This 
is just the standard BTK contribution [2, 34]. The non-trivial contribution is 
represented by the last term in Eq. (1.59) which accounts for the presence of 
the second NS interface. We observe that this non-local contribution to the 
current is small as oc IppL^. This term will be analyzed in details below. 

The functions and 5^„t are the Keldysh Green functions evaluated 
on the trajectories 1*" and 1°"* respectively. Using the above expression for 
the Riccati amplitudes and the distribution functions we find 

Sp(f3ff^$„* - f3ft^„) = 2 Spir^ir^i" - f/'f/'^ - l](/io(£ + eVi) - hois)) 
- 2 Sp[f^ir^; - - l]{ho{e - eV^) - ho{s)) 

+ 2(1 - tanh^ iLQ/vp) Sp[K^^{ST2Si,^{hQ{e + eV2) - ho{e)) 
+ \a\''S2'2'(T2S+^,S^2'<^2S+^,{ho{e - eV2) - ho{e))}K+~' 

X {S+,Siy - \a''\S+,,a2SyiS+,a2Syy)], (1.61) 

where we explicitly used the fact that in equilibrium Sp(f35^„t — rsg^in) = 0. 
Substituting (1.61) into (1.59), we finally obtain 

h = ir^'iVi) + hi{Vi) + Ii2{V2). (1.62) 

The correction to the local BTK current (arising from trajectories crossing 
also the second NS interface) has the following form 

hiiVi) = ~^J de{SpKr^^ - f /^A''^ - mo{s + eV,) - ho{s)) 

- Mr?ir?i" - rfTf + - l]{ho{e - eV,) - ho{e))}, (1.63) 
while for the cross-current we obtain 
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112(^2) = J de{l - tanh^ iLH/vp) 

X Sp[K^^^{S2'2S+^{ho{e + eV2) - /10(e)) 
+ \a\''S2'2"J2S+^,S^2''^2S+^,ihoie - e^a) - hoie))}K+~" 

X {S+,Sn' - \a''\S+,,a2Si,iSt,,a2Srv)]. (1.64) 

Eqs. (1.62)-(1.64) fully determine the current across the first interface at 




Fig. 1.5. Diagrams representing four different contributions to the cross-current 
I12 (1.64). Solid (dotted) lines correspond to propagating electron-like (hole-like) 
excitations and t2i = / cosh.{iLQ /vf)- 

arbitrary voltages, temperature and spin-dependent interface transmissions. 

In right hand side of Eq. (1.64) we can distinguish four contributions with 
different products of S-matrices. Each of these terms corresponds to a cer- 
tain sequence of elementary events, such as transmission, reflection, Andrcev 
reflection and propagation between interfaces. Diagrammatic representation 
of these four terms is offered in Fig. 1.5. The amplitude of each of the pro- 
cesses is given by the product of the amplitudes of the corresponding ele- 
mentary events. For instance, the amplitude of the process in Fig. 1.5c is 
/ = — iS'ii/t2iS'2'2'(icr2^22'- Eq.(1.64) this process is identified by the term 
Sp(//"'") with the hole distribution function as a prefactor. It is straightfor- 
ward to observe that the processes of Fig. 1.5a, 1.5b and 1.5d correspond to 
the other three terms in (1.64). We also note that the processes of Fig. 1.5a and 
1.5d describe direct electron (hole) transport, while the processes of Fig. 1.5b 
and 1.5c correspond to the contribution of CAR. 

Assuming that both interfaces possess inversion symmetry as well as re- 
flection symmetry in the plane normal to the corresponding interface we can 
choose iS-matrices in the following form 
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Sii = Si'v = = Sj^ii, = U{ip) 







/2 



-^0^/2 I U+{<^), 



(1.65) 



S22 = S2'2' = S_22 = S_2'2' ~ 
S22' = S2'2 = S_22i = S_2i2 = i 



V:R^e-«=^/2 

/D^e^^=/2 

/D^e-^^=/2 



(1.66) 
(1.67) 

(1.68) 



Here -Ri(2)t(4,) = 1 ~ -Di(2)t(^) are the spin dependent reflection coefflcients 
of both NS interfaces, 6\^2 are spin-mixing angles and U{(f) is the rotation 
matrix in the spin space which depends on the angle </? between polarizations 
of the two interfaces, 



U{^p) = exp{—i(pai/2) 



cos(ip/2) — isin((/5/2)\ 
— isin((^/2) cos(<^/2) J ' 



(1.69) 



In general spin current is not conserved in heterostructures with spin active 
interfaces. However single barrier with S'-matrix (1.67)- (1.68) does not violate 
spin current conservation [40]. It is easy to show that in our two barrier 
structure with interface 5-matrices (1.65)-(1.68) spin current is conserved in 
the normal state for arbitrary barriers polarizations and in superconducting 
state for collinear barriers polarizations. 

Substituting the above expressions for the ^-matrices into Eqs. (1.63) and 
(1.64) we arrive at the flnal results for both /ii(Fi) and /i2(V2) which will be 
specified further below. 



1.2.6 Cross-current 

First let us consider the cross-current /i2(V2). From the above analysis we 
obtain 



h2{V2] 



Go 
' 4e 



/ 



ds 



tanh ■ 



tanh ■ 



1 — tanh^ iLfi/vp 
W{zi,Z2,e,ip) 



2T 2T 

X { [Di^D2i - |a|2Di^D2;(J?i^ + J?2t) + lal^-Du-Rit-Da^-Rat] \K{zu ^2, £)|'c 

+ [^it^2t - |a|2£)^£)2t(^it + R2i) + |a|4Diti?uI>2t^2j \K{zl, z*2,e)\^c 
+ [Di^D2i - \afDnD2i{Rit+R2t) + \a\^Di^Ri^D2iR2t] \K{zl,Z2,e)\''s 

+ [DiiD2t - |a|2Ditr>2t(i?u + ^2;) + \a\^DnRi^D2^R2i\ \K{zuzle)\''sY 

(1.70) 
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where we define c = cos^((/?/2), s = sin^((/3/2), 

K{zi,Z2, £) = (!- 0^2:122) - [e(l + 0^2122) + ^a(2i + Z2)] Q, (1.71) 
W{zi,Z2,e,^) = \K{zi,Z2,e)K{zl,z;,e)cos'^{ip/2) 

+Kizlz2,e)Kizi,z;, s) sm'{if/2)\ 



2. /OM2 (1-72) 



and Zi = -^/%^exp(i6li) ( i = 1, 2). 

Eq. (1.70) represents our central result. It fully determines the non-local 
spin-dependent current in our three-terminal ballistic NSN structure at arbi- 
trary voltages, temperature, interface transmissions and polarizations. 

Let us introduce the non-local differential conductance 

Before specifying this quantity further it is important to observe that in gen- 
eral the conductance Gi2(V^) is not an even function of the applied voltage ¥2- 
This asymmetry arises due to formation of Andreev bound states in the vicin- 
ity of a spin-active interface [41, 42]. It disappears provided the spin mixing 
angles 9i and O2 remain equal to or tt. 

In the normal state we have /i2(V2) = —GN12V2, where 

Gn,, = ^ [{DuD2i + Di^D2t) cos\ip/2) 

+ {Di^D2i + Di^D2t) sin2(<^/2)] . (1.74) 

Turning to the superconducting state, let us consider the limit of low 

temperatures and voltage T, V2 <C A. In this limit only subgap quasiparti- 
cles contribute to the cross-current and the differential conductance becomes 
voltage-independent, i.e. I12 = —G12V2, where 



G12 = Go(l - tanh^ LA/vf) 



Di^DiiD2^D2i 



[\K{zi,Z2,0Wd+\K{zi,z^,0)\^s 

^ ^'(|i^(zi,Z2,0)P5+|if(zi,z*,0)p,?)2J 

where, as before, c — cos^((/^/2) and s = sin^(iyj/2). In the case of spin-isotropic 
interfaces Eqs. (1.75) and (1.70) reduce to the results [14]. 

Provided at least one of the interfaces is spin-isotropic. the conductance 
(1.75) is proportional to the product of all four transmissions Di'^Dx\^D2\D2\., 
i.e. it differs from zero only due to processes involving scattering with both spin 
projections at both NS interfaces. As in the case of spin-isotropic interfaces 
[14] the value G12 (1.75) gets strongly suppressed with increasing L, and at 
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Fig. 1.6. Zero temperature differential non-local conductance as a function of volt- 
age at zero spin-mixing angles 61^2 = 0. 

sufRciently high interface transmissions this dependence is in general non- 
exponential in L. In the spin-degenerate case for a given L the non-local 
conductance reaches its maximum for reflectionless barriers D12 — 1- In this 
case we arrive at a simple formula 

G12 = Go{l-t£inh^ LA/vf). (1.76) 

We observe that for small L <C vp/A the conductance G12 identically coin- 
cides with its normal state value GjVia = Gq at any temperature and voltage 
[14]. This result implies that CAR vanishes for fully open barriers. Actually 
this conclusion is general and applies not only for small but for any value of 
L, i.e. the result (1.76) is determined solely by the process of direct electron 
transfer between N-terminals for all L. 

At the first sight, this result might appear counterintuitive since the be- 
havior of ordinary (local) AR is just the opposite: It reaches its maximum 
at full barrier transmissions. The physics behind vanishing of CAR for per- 
fectly trasparent NS interfaces is simple. One observes (cf. Fig. 1.1) that CAR 
inevitably implies the flow of Cooper pairs out of the contact area into the su- 
perconducting terminal. This flow is described by electron trajectories which 
end deep in the superconductor. On the other hand, it is obvious that CAR 
requires "mixing" of these trajectories with those going straight between two 
normal terminals. Provided there exists no normal electron reflection at both 
NS interfaces such mixing does not occur, CAR vanishes and the only re- 
maining contribution to the non-local conductance is one from direct electron 
transfer between N-terminals. 
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This situation is illustrated by the diagrams in Fig. 1.5. It is obvious that in 
the case of non-reflecting NS interfaces only the process of Fig. 1.5a survives, 
whereas all other processes (Fig. 1.5b, 1.5c and 1.5d) vanish for reflectionless 
barriers with i?i(2)t(4.) = 0. The situation changes provided at least one of the 
transmissions is smaller than one. In this case scattering at SN interfaces mixes 
up trajectories connecting Ni and N2 terminals with ones going deep into and 
coming from the superconductor. As a result, all four processes depicted in 
Fig. 1.5 contribute to the cross-current and CAR contribution to 6*12 does 
not vanish. 

In the limit |eV2|, T <C A and at zero spin-mixing angles 61^2 — from Eq. 
(1.75) we obtain 

- |if(zi,Z2,0)P l^it^U^2T^2i 

+ (At-Z?u)(C2T-i52;)cos(p}. (1.77) 

In the lowest (first order) order in the transmissions of both interfaces and 
for collinear interface polarizations Eq. (1.77) reduces to the result by Falci 
et al. [13] provided we identify the tunneling density of states NgDif, NqDi^, 
A*oi?2ti ^-nd NqD2i with the corresponding spin-resolved densities of states 
in the ferromagnetic electrodes. For zero spin-mixing angles and low voltages 
the L-dependence of the nonlocal conductance G12 reduces to the exponential 
form G12 c>c exp{~2LA/vp) either in the limit of small transmissions or large 
L > vf/A. 
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At arbitrary voltages and temperatures the cross-current has a simple 
dependence in the limit of zero spin mixing angles {61^2 = 0) 



h2{v, V2) = Ii2{v = 0, V2) cos2(<p/2) + h2{v = TT, V2) sm\ip/2), (1.78) 



i.e. in this limit at any ip the nonlocal current is equal to a proper superposi- 
tion of the two contributions corresponding to parallel (i^ = 0) and antipar- 
allel {(fi = it) interface polarizations. Some typical curves for the differential 
non-local conductance are presented in Fig. 1.6 at sufficiently high interface 
transmissions and zero spin mixing angles 61^2 = 0. 

Let us now turn to the limit of highly polarized interfaces which is ac- 
counted for by taking the limit of vanishing spin-up (or spin-down) transmis- 
sion of each interface. In this limit our model describes an HSH structure, 
where H stands for fully spin-polarized half-metallic electrodes. In this case 
we obtain (£>i-|- = -Di, -D14, = 0, D2f = D2, and D21 = 0) 



We observe that the nonlocal conductance has opposite signs for parallel {(p = 
0) and antiparallel {(p = tt) interface polarizations. We also emphasize that, as 
it is also clear from Eq. (1.77), the cross-conductance G12 of HSH structures - 
in contrast to that for NSN structures - does not vanish already in the lowest 
order in barrier transmissions -Di-|--D2t- 

In general the non-local conductance is very sensitive to particular values 
of the spin-mixing angles 61 and 62, as illustrated, e.g., in Fig. 1.7. Comparing 
the voltage dependencies of the nonlocal conductance evaluated for the same 
transmissions and presented in Figs. 1.6 and 1.7, we observe that they can 
differ drastically at zero and non-zero values of 9i^2- 

At low voltages and temperatures and at zero spin mixing angles the non- 
local conductance of HSH structures is determined by Eq. (1.77) with Di^ = 
D21 = 0. For fully open barriers (for "spin-up" electrons) Di^ = D2f = 1 we 
obtain 



Interestingly, for = this expression exactly coincides with that for fully 
open NSN structures, Eq. (1.76). At the same time for small L the result 
(1.80) turns out to be 2 times bigger that the analogous expression in the 
normal case, i.e. for fully open HNH structures, cf. Eq. (1.74). This result 
can easily be interpreted in terms of diagrams in Fig. 1.5. We observe that - 
exactly as for the spin degenerate case - CAR diagrams of Fig. 1.5b,c vanish 
for reflectionless barriers, whereas diagrams of Fig. 1.5a,d describing direct 
electron transfer survive and both contribute to G12. Thus, CAR vanishes 
identically also for fully open HSH structures. The factor of 2 difference with 
the normal case is due to the fact that the diagram of Fig. 1.5d vanishes in 
the normal limit. 




X 



{[l + \af] \Kizlz*2,e)\^c-2\a\^\Kiz^,z*2,e)\^~s}. (1.79) 



G12 = Go(l - tanh^ LA/vp) cosip. 



(1.80) 
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1.2.7 Correction to BTK 

Using the above formalism one can easily generalize the BTK result to the 
case of spin-polarized interfaces [34]. For the first interface we have 

ir^'iVi) = ^J + eFi) - ho{e)]{l + |a|2) 

Here transmission and reflection coefficients as well as the spin mixing angle 
depend on the direction of the Fermi momentum. In the spin-degenerate case 
the above expression reduces to the standard BTK result [2] . 

Evaluating the nonlocal correction to the BTK current due to the presence 
of the second interface we arrive at a somewhat lengthy general expression 

hi{Vi) = ^ [ de{ho{e + eV^)-ho{e))-—^ -{2W{zi, Z2,e,^) 

ze J w [zi,z2, e, f) 

- Ri^\cK{zi/Ri-f,z2,e)K{zl,z2,e) + sK{zi/Rif,z2,e)K{zl,z2,e)\^ 

- Rii\cK{zl/Rii, z2,e)K{zi,Z2, e) + sK{zl/Rii, Z2, e)K{zi,z2.,e)\^^ 

+ ^ [deiho{e + eV,)-ho{e))-^^l^{ 

4e J W{zi,Z2,s,ip) L 

|a|^|cif (0, Z2, s)K{zl, Z2,e) + sK{0, z'2,e)K{zl, 22, e)|^ 
+ |ap|cX(0, z;,e)K{zi,Z2, s) + sK{0, Z2, e)Kizi, Z2,s)\^ 

+ -^\cK'{z;,e)K{zi,Z2,e) + sK'{z2,e)K{zi,z;,e)\'' 
\a\ 

+ j^\cK'{z2,e)K{zlz*2,s) + sK'{z*2,s)K{zlz2,e)\^} 

+ ^R2^R2ism\e2/2)Sc 



j de{ho{e + eVi) - ho{e)) 



(1 - tanh^ iLfl/vpf 



W{zi,Z2,s,ip) 

X [\af{Dl^ + Dl) - 2|a|^£>it£>u(i?it + i?u) + \af{Dl^Ru + D^RI^)] , 

(1.82) 

where K'{z2,e) = dK{zi,Z2,£)/dzi. This expression gets significantly simpli- 
fied in the limit of zero spin- mixing angles ^1,2 = in which case we obtain 
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/ii(^i) = § / deihoie + eV^) - 



(£)) 



""^^ |i^(zi,.2,s)P '''' \K{z,,z,,eW 

mz,,z,,ew /• ^'-''^ 

In contrast to the expression for the cross-current I12 (cf. Eq. (1.78)), in the 
limit of zero spin- mixing angles the correction In to the BTK current does 
not depend on the angle between the interface polarizations. In particular, 
at |eVi|,T <C /\ we have In = GnVi where 



[1 + Z1Z2 + {zi + Z2) tanhLA/vp]'^ 
(1 + Z2 tav^LA/vpf + 3(^2 + tar^LA/vpf 
+ [l + z^Z2 + {z,+Z2)te.nhLA/vp]^ • ^^'^^^ 

In the tunneling limit Dif, Di^, D2f, -D24, <C 1 we reproduce the result of Ref. 
[13] 

Gn = ^(I?it + D,^)iD2^ + D21) eM-^LA/vp), (1.85) 

which turns out to hold at any value (p. 

As compared to the BTK conductance the CAR correction (1.82) contains 
an extra small factor A2/L'^ and, hence, in many cases remains small and can 
be neglected. On the other hand, since CAR involves tunneling of one electron 
through each interface, for strongly asymmetric structures with Dif,Dii <C 
1 and D2i;,D2_i ^ 1 it can actually strongly exceed the BTK conductance. 
Indeed, for Di^^ <C 1, -R2t-^24, 1 and provided the spin mixing angle 9i is 
not very close to tt from Eq. (1.82) we get 

n = G,ADi^ + DK) .^8g^ 

cosh(2Lzi/wir) + cos6'isinh(2I,Z\/wir)' ^ ' ' 



i.e. for 



' — < -r2-exp(-2LZ\/t;F) 



the contribution (1. 86) may well exceed the BTK term Gf^^ cx Dif Di^. The 
existence of such a non-trivial regime further emphasizes the importance of 
the mechanism of non-local Andreev reflection in multi-terminal hybrid NSN 
structures. 



1.3 Diffusive FSF structures 



Lot us now turn to the effect of disorder. In what follows we will consider a 
three-terminal diffusive FSF structure schematically shown in Fig. 1.8. Two 
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ferromagnetic terminals Fi and F2 with resistances r^Vi and r^Vs and electric 
potentials Vi and V2 are connected to a superconducting electrode of length 
L with normal state (Drude) resistance and electric potential ^ = via 
tunnel barriers. The magnitude of the exchange field hi^2 = |^i,2| in both fer- 
romagncts Fi and F2 is assumed to be much bigger than the superconducting 
order parameter A of the S-terminal and, on the other hand, much smaller 
that the Fermi energy, i.e. A <^ hi^2 ep- The latter condition allows to 
perform the analysis of our FSF system within the quasiclassical formalism of 
Usadel equations for the Green- Keldysh matrix functions G formulated below. 




Fig. 1.8. FSF structure under consideration. 



1.3.1 Quasiclassical equations 

In each of our metallic terminals the Usadel equations can be written in the 
form [32] 

iDV(GVG) = [f2 + eV, G], & = 1, (1.87) 

where D is the diffusion constant, V is the electric potential, G and fi are 
8x8 matrices in Keldysh-Nambu-spin space (denoted by check symbol) 

= (1.89) 

£ is the quasiparticle energy, A{T) is the superconducting order parameter 
which will be considered real in a superconductor and zero in both fcrromag- 
nets, h = hi(2) in the first (second) ferromagnetic terminal, h = outside 
these terminals and & = [ai, 0'2, cts) are Pauli matrices in spin space. 

Retarded and advanced Green functions G^ and G^ have the following 
matrix structure 

G^'^ = [%,A ^^r,a) ■ (1-90) 
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Here and below 2x2 matrices in spin space are denoted by hat symbol. 

Having obtained the expressions for the Grcen-Keldysh functions G one 
can easily evaluate the current density j in our system with the aid of the 
standard relation 

3 = j Sp[r3(GVG)^]rf£, (1.91) 

where a is the Drude conductivity of the corresponding metal and rs is the 
Pauli matrix in Nambu space. 

In what follows it will be convenient for us to employ the so-called Larkin- 
Ovchinnikov parameterization of the Keldysh Green function 

G^ = G«/-/G^, f = fL + r3fT, (1.92) 

where the distribution functions Jl and /t are 2x2 matrices in the spin 
space. 

For the sake of simplicity we will assume that magnetizations of both 
ferromagnets and the interfaces (see below) are coUinear. Within this approx- 
imation the Green functions and the matrix f2 are diagonal in the spin space 
and the diffusion-like equations for the distribution function matrices /l and 
/t take the form 

-£>V (£'^(r, e)V/T(r, s)) + 2E{r, e)fT{r, s) = 0, (1.93) 

-£»V (l)^(r, e)V/L(r, s)) = 0, (1.94) 

where 

S{r, e) = -lA Im F^, (1.95) 

= (Re G^) ' + (im P") ' , (1.96) 

= (ReG^)^ - (ReF'^y . (1.97) 

The function S{r,e) differs from zero only inside the superconductor. It ac- 
counts both for energy relaxation of quasiparticles and for their conversion to 
Cooper pairs due to Andreev reflection. The functions and acquire 
space and energy dependencies due to the presence of the superconducting 
wire and renormalize the diffusion coefficient D. 

The solution of Eqs. (1.93)-(1.94) can be expressed in terms of the diffuson- 
like functions and V'" which obey the following equations 



-DV 



£)^(r,£)VP^(r,r',£) +2IJ{r,e)V^{r,r',e) = S{r - r'), (1.98) 



L»^(r,e)VP^(r,r',£) = 6{r - r'). (1.99) 
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1.3.2 Boundary conditions 

The solutions of Usadel equation (1.87) in each of the metals should be 
matched at FS-interfaces by means of appropriate boundary conditions which 
account for electron tunneling between these terminals. The form of these 
boundary conditions essentially depends on the adopted model describing 
electron scattering at FS-interfaces. As before, we stick to the model of the 
so-called spin-active interfaces which takes into account possibly different bar- 
rier transmissions for spin-up and spin-down electrons. Here we employ this 
model in the case of diffusive electrodes and also restrict our analysis to the 
case of tunnel barriers with channel transmissions much smaller than one. In 
this case the corresponding boundary conditions read [43, 44] 

Aa+G+d,G+ = ^[G_,G+] + ^[{amTs,G_},G+] + i^[&mT3,G+], 

(1.100) 

-^<7_G_a,G_ = ^[G+,G-] + ^[{<Tmr3,G+},G-] +z^[<Tmr3,G-], 

(1.101) 

where G- and G+ are the Grecn-Kcldysh functions from the left (x < 0) 
and from the right {x > 0) side of the interface, A is the effective contact 
area, m is the unit vector in the direction of the interface magnetization, 
a± are Drude conductivities of the left and right terminals and Gt is the 
spin-independent part of the interface condiictancc. Along with Gt there also 
exists the spin-sensitive contribution to the interface conductance which is 
accounted for by the G^-term. The value Gm equals to the difference between 
interface conductances for spin-up and spin-down conduction bands in the 
normal state. The Gi^-term arises due to different phase shifts acquired by 
scattered quasiparticles with opposite spin directions. 

Employing the above boundary conditions we can establish the following 
linear relations between the distribution functions at both sides of the interface 

Aa+DldJ+T = Aa-D^dJ-T = 5t(/+t - f-r) + 9m{f+L - I-l), 

(1.102) 

Acj+b%dJ+L = Aa_b^_dJ_L = hif+L - f-L)+9m{f+T ' /-t), 

(1.103) 

where Qt, Ql, and Qm are matrix interface conductances which depend on the 
retarded and advanced Green functions at the interface 

Qt = Gt (RcG^) (RcG^) + (imFf ) (imF^^) , (1.104) 

= Gt [{ReG+) (RcG^) - (RcF^) (RbFj^)] , (1.105) 

= GmO-m (Re G^) (rc G^!) . (1.106) 
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Note that the above boundary conditions for the distribution functions do not 
contain the Gi^-tcrm cxpHcitly since this term in Eqs. (I.IOO)-(I.IOI) does not 
mix Green functions from both sides of the interface. 

The current density (1.91) can then be expressed in terms of the distribu- 
tion function /t as 

i = - f / Sp[b^VfT]ds. (1.107) 



1.3.3 Spectral conductances 

Let us now employ the above formahsm in order to cvakiatc electric currents 
in our FSF device. The current across the first (SFi) interface can be written 
as 

/i = ^ y giiie) [fo{£ + eV,) - /o(e)] de 

- \ j 9i2{£) [h{e + eV2) - Ms)] de, (1.108) 

where /o(e) = tanh(e/2r), gu and 512 are local and nonlocal spectral electric 
conductances. Expression for the current across the second interface can be 
obtained from the above equation by interchanging the indices 1^2. Solving 
Eqs. (1.93)-(1.94) with boundary conditions (1.102)-(1.103) we express both 
local and nonlocal conductances guj (s) in terms of the interface conductances 
and the function V. The corresponding results read 

gui^) = {R2-^^ + ^2^2^1m — -Rf-Rim + ■Rl2-^12-^2m — RlmR^m)^! 

(1.109) 

512(e) = 521(e) = (.Rfj-^^ + -^2 -^12-^lm + -Rl2-^lm-R2m + -Ri2-^l'-^2m)^, 

(1.110) 

where we defined 

M^''^ = kV'kr-{Rl,')\ (1.111) 

K,~^ = M^M^ + Rlm^lm ~ ^2^2^1771 ^ 2i?f2-Kl2-^lm-R2m - Rl Rl Rim 

(1.112) 

and introduced the auxiliary resistance matrix 

R( = 5iT(e)[5iT(e)5iL(e) - 'gl^{£)]~^ 

(Ti as 

The resistance matrices R2, R^ and R2 can be obtained by interchanging the 
indices 1 o 2 and T ^ L in Eq. (1.113). The remaining resistance matrices 
Rj^2 ^'^d Rjm are defined as 
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^T^L = ^T,. ^ DsVl'\rur2,e) ^ 

Rjm = gjm{e)[9jT{e)gjL{£) - 5|m(£)]"\ (1.115) 

where j = 1,2. The spectral conductance gij can be recovered from the matrix 
gij simply by summing up over the spin states 

9iM = lsp[gij{e)]. (1.116) 

It is worth pointing out that Eqs. (1.109), (1.110) defining respectively 
local and nonlocal spectral conductances are presented with excess accuracy. 
This is because the boundary conditions (I.IOO)-(I.IOI) employed here remain 
applicable only in the tunneling limit and for weak spin dependent scattering 
<C! Gt- Hence, strictly speaking only the lowest order terms in 
2 and Gy,^ ^ need to be kept in our final results. 
In order to proceed it is necessary to evaluate the interface conductances 
as well as the matrix functions '£^i'2s- Restricting ourselves to the second 
order in the interface transmissions we obtain 

gMs) = Gr^Osirue) + G%,^^^^H^), (1.117) 

giL{e) = GT,%(ri,£) - G|^^^Jr^^t>i(e), (1-118) 
gim{£) = Gmii's{ri,e)&mi, (1.119) 

and analogous expressions for the interface conductances of the second inter- 
face. The matrix function 



Ui (e) = ^ {Re [Ci (n , ri , 2/1+) + Ci (n , n , 2h^ )] 



- o-mi Re [Ci(ri, n, 2/1+) - Ci(ri, n, 2/1^)] } (1.120) 

with hf = /ii ± £ defines the correction due to the proximity effect in the 
normal metal. 

Taking into account the first order corrections in the interface transmis- 
sions one can derive the density of states inside the superconductor in the 
following form 



me^-A^) ^ Ds ^ 
vs{r,e) = , =^ H p5 5- > 



GT,ReCs(r,ri,2w^) 



- &miG,p^ ImCs(r, r^, 2u)^) 



(1.121) 



where 
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Ve^ - e> A, 

oj^ = { iVA^-e^, \e\ < A, (1.122) 
^-Vs^-A^, s<A, 

and the Cooperon Cj{r,r',s) represents the solution of the equation 

(-DV^ - ie) C{r, r', e) = 5{r - r') (1.123) 

in the normal metal leads (j = 1, 2) and the superconductor {j = 5). In the 
quasi-one-dimensional geometry the corresponding solutions take the form 

C,(x„x„e) = '''f^^^''^\ J = 1,2, (1.124) 

, sinh[fcs(L-x')]sinhfcsa; , /i iot;\ 

Cs{x,x,e) = — , „ „ ■ ufi. r\ — ' a; >a;, 1.125 



where <S's,i,2 are the wire cross sections and fci^2,s ~ ie/-Di,2,s- 

Substituting Eq. (1.121) into Eqs. (1.117) and (1.118) and comparing the 
terms a G^^ we observe that the tunneling correction to the density of states 
dominates over the terms proportional to Ui which contain an extra small 
factor sj Ajh <C 1. Hence, the latter terms in Eqs. (1.117) and (1.118) can be 
safely neglected. In addition, in Eq. (1.121) we also neglect small tunneling 
corrections to the superconducting density of states at energies exceeding the 
superconducting gap A. Within this approximation the density of states inside 
the superconducting wire becomes spin-independent !>5(r,£) = cFQVs{r,E). It 
can then be written as 

+ GT.ReC5(r,r.,2a,''). (1.126) 

Accordingly, the interface conductances take the form 

5iT(e) = giL{e) = GT,vs{ri,e), (1.127) 
5im(e) = Grmi's{ri,e)ami. (1.128) 

Let us emphasize again that within our approximation the G<p-term does 
not enter into expressions for the interface conductances (1.127)-(1.128) and, 
hence, does not appear in the final expressions for the conductances gij{s). 

In the limit of strong exchange fields hi ^2 ^ A and small interface 
transmissions considered here the proximity effect in the ferromagnets re- 
mains weak and can be neglected. Hence, the functions P^'^(ri, ri, e) and 
■D^'^(r2,r2,e) can be approximated by their normal state values 
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Pf'^(ri,ri,e) = airjv.iM, (1.129) 

't>l^''{r2,r2,e)=(J2rN^/D2, (1.130) 
rN,=h/{ajSj), j = l,2, (1.131) 

where rjvi and rjva are the normal state resistances of ferromagnetic termi- 
nals. In the the superconducting region an effective expansion parameter is 
Gti_2^Cs(^)' where r^g{e) = ^s{£)/{<^sSs) is the Drude resistance of the su- 
perconducting wire segment of length ^s:(e) = ■\/Ds/2\uj^\ and is the 
function of e according to Eq. (1.122). In the limit 

GT,,.r^,(e)«l, (1.132) 

which is typically well satisfied for realistic system parameters, it suffices to 
evaluate the function 'Dg{x,x' ,e) for impenetrable interfaces. In this case we 
find 



V^{x,x',e)={ 



-Csix,x',0), \e\>A. 



Wc note that special care should be taken while calculating Pg(a;, x',£) at 
subgap energies, since the coefficient in Eq. (1.94) tends to zero deep inside 
the superconductor. Accordingly, the function 'Dg{x, x' , e) becomes singular in 
this case. Nevertheless, the combinations Rj'{A4^)~^ and Ri2{M.^)~^ remain 
finite also in this limit. At subgap energies we obtain 



1 



2Ke<^/€s(^) 



(1.134) 



where k = 1 — e^/Z\^ and d = \x2 — Xi\ is the distance between two FS 
contacts. Substituting the above relations into Eq. (1.110) we arrive at the 
final result for the non-local spectral conductance of our device at subgap 
energies (|e| < A) 

„ i^nsi^)exp[-d/^s{e)] 

9i2{£) — g2i{£) — 



1 + 



2[rjvi + l/5Ti(£)][riv. + 1/5t2(£)] 

mim2 Gml Gm2 1 

K 9Ti{e) gT2{e) ^ + + '^^V ;,(£)GT,GT.e-^/^^-(^) 

(i".135) 

Eq. (1.135) represents the central result of this section. It consists of two 
different contributions. The first of them is independent of the interface po- 
larizations mi^2. This term represents direct generalization of the result [24] 
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in two different aspects. Firstly, the analysis [24] was carried out under the 
assumption r^Vi 2.9ti,2(£) ^ 1 which is abandoned here. Secondly (and more 
importantly), sufficiently large exchange fields /ii_2 ^ ^ of ferromagnetic elec- 
trodes suppress disorder-induced electron interference in these electrodes and, 
hence, eliminate the corresponding zero-bias anomaly both in local [27, 28, 29] 
and non-local [24] spectral conductances. In this case with sufficient accuracy 
one can set gni^) = GTii^s{xi,e) implying that at subgap energies gnis) is 
entirely determined by the second term in Eq. (1.126) which yields in the case 
of quasi-one-dimensional electrodes 



gT2{e) = 



2(Z\2-e2) 

2(zi2 _ £2) 



Gt^ + Gr.e-'^^^'''-^'^ 



(1.136) 
(1.137) 



Note, that if the exchange field /ii^2 in both normal electrodes is reduced 
well below A and eventually is set equal to zero, the term containing Ui{e) in 
Eqs. (1.117), (1.118) becomes important and should be taken into account. In 
this case we again recover the zero-bias anomaly [27, 28, 29] gni^) oc 
and from the first term in Eq. (1.135) we reproduce the results [24] derived in 
the limit hi^2 — >■ 0. 

The second term in (1.135) is proportional to the product m\m2GmiGm2 
and describes non-local magnctoconductancc effect in our system emerging 
due to spin-sensitive electron scattering at FS interfaces. It is important that 
- despite the strong inequality \Gmi\ «C Gn - both terms in Eq. (1.135) can 
be of the same order, i.e. the second (magnetic) contribution can significantly 
modify the non-local conductance of our device. 

In the limit of large interface resistances rjVi,2fl'Ti,2(£) 1 the formula 
(1.135) reduces to a much simpler one 



512(e) = 521(e) = TiM e^^[-d/^s{e)\ 

\2 - £2 A 

gTl{e)gT2{e) + mim2GmlGm2 ^2 _ g2 



(1.138) 



Interestingly, Eq. (1.138) remains applicable for arbitrary values of the angle 
between interface polarizations m\ and m2 and strongly resembles the anal- 
ogous result for the non-local conductance in ballistic FSF systems (cf., e.g., 
Eq. (1.77) in the previous section). The first term in the square brackets in 
Eq. (1.138) describes the fourth order contribution in the interface transmis- 
sions which remains nonzero also in the limit of the nonfcrromagnctic leads 
[24]. In contrast, the second term is proportional to the product of transmis- 
sions of both interfaces, i.e. only to the second order in barrier transmissions. 
This term vanishes identically provided at least one of the interfaces is spin- 
isotropic. 
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Fig. 1.9. Local (long-dashed line) and non-local (short-dashed and solid lines) spec- 
tral conductances normalized to its normal state values. Here we choose vni = rjvj = 
5r^^{0), xi = L - X2 ^ 5^5 (0), X2 ~ xi ^ Cs(0), Gti = Gt2 = 4Gmi = 46^2 = 
0.2/r^g (0). Energy dependence of non-local conductance is displayed for parallel (P) 
mim2 = 1 and antiparallel (AP) mim2 = —1 interface magnetizations. Inset: The 
same in the limit of low energies. 



Contrary to the non-local conductance at subgap energies, both local con- 
ductance (at all energies) and non-local spectral conductance at energies above 
the superconducting gap are only weakly affected by magnetic effects. Neglect- 
ing small corrections due to term in the boundary conditions we obtain 

gn{e)^R^{M^r\ 922{e) ^ {M^)-\ (1.139) 
512(e) = .921(e) = Rf^iM^y^ \e\ > A. (1.140) 

Eqs. (1.139) and (1.140) together with the above expressions for the non- 
local subgap conductance enable one to recover both local and non-local spec- 
tral conductances of our system at all energies. Typical energy dependencies 
for both 5ii(e) and 512(e) are displayed in Fig. 1.9. For instance, we observe 
that at subgap energies the non-local conductance 1712 changes its sign being 
positive for parallel and negative for antiparallel interface polarizations. 

1.3.4 I-V curves 

Having established the spectral conductance matrix gij{e) one can easily re- 
cover the complete I —V curves for our hybrid FSF structure. In the limit of 




-0.2 I ' ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 

Fig. 1.10. Non-local resistance (normalized to its normal state value) versus tem- 
perature (normalized to the superconducting critical temperature Tc) for parallel 
(P) and antiparallel (AP) interface magnetizations. The parameters are the same as 
in Fig. 1.9. 




- 



-0.5 I ' ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 

Fig. 1.11. The same as in Fig. 1.10 for the following parameter values: vmi = ~ 

5r53(0), x^=L-x2 = 5Cs(0), X2 - = Cs(0), = Gt^ = 25G„i = 250^^ = 
0.025/res(0). 
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low bias voltages these I — V characteristics become linear, i.e. 

h = Gu{T)Vi - Gi2{T)V2, (1.141) 
h = -G2i{T)Vi + G22{T)V2, (1.142) 

where G^- (T) represent the linear conductance matrix defined as 

GijiT) = ^ fgijie)^^. (1.143) 
4J J cosh^ — 

It may also be convenient to invert the relations (1.141)-(1.142) thus express- 
ing induced voltages Vi^2 in terms of injected currents Ii^2- 

Vi = Rii{T)h + Ri2(T)l2, (1.144) 

V2 = R2l{T)h + i?22(T)/2, (1.145) 

where the coefficients Rij{T) define local {i = j) and nonlocal {i ^ j) resis- 
tances 

^"'^'° Gn(T)G°m-gyr) ' 

^"'^'°^"(^'° Gn(T)G°'m-Gj,(T) 

and similarly for i?22(r). In non- ferromagnetic NSN structures the low tem- 
perature non-local resistance R\2 {T — >■ 0) turns out to be independent of both 
the interface conductances and the parameters of the normal leads [24] . How- 
ever, this universality of if!i2 does not hold anymore provided non-magnetic 
normal metal leads are substituted by ferromagncts. Non-local linear resis- 
tance Ri2 of our FSF structure is displayed in Figs. 1.10, 1.11 as a function 
of temperature for parallel {mim2 = 1) and antiparallel {m,\vn,2 = — 1) in- 
terface magnetizations. In Fig. 1.10 we show typical temperature behavior of 
the non-local resistance for sufficiently transparent interfaces. For both mu- 
tual interface magnetizations R\2 first decreases with temperature below Tc 
similarly to the non-magnetic case. However, at lower T important differences 
occur: While in the case of parallel magnetizations R\2 always remains pos- 
itive and even shows a noticeable upturn at sufficiently low T, the non-local 
resistance for antiparallel magnetizations keeps monotonously decreasing with 
T and may become negative in the low temperature limit. In the limit of very 
low interface transmissions the temperature dependence of the non-local resis- 
tance exhibits a well pronounced charge imbalance peak (see Fig. 1.11) which 
physics is similar to that analyzed in the case of non-ferromagnetic NSN struc- 
tures [18, 22, 24]. Let us point out that the above behavior of the non-local 
resistance is qualitatively consistent with available experimental observations 
[5] . More experiments would be desirable in order to quantitatively verify our 
theoretical predictions. 
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1.4 Concluding remarks 

In this paper we developed a non-perturbative theory of non-local electron 
transport in both ballistic and diffusive NSN and FSF three-terminal struc- 
tures with spin-active interfaces. Our theory is based on the quasiclassical 
formalism of energy-integrated Green-Eilenberger functions supplemented by 
appropriate boundary conditions describing spin-dependent scattering at NS 
and FS interfaces. Our approach applies at arbitrary interface transmissions 
and allows to fully describe non-trivial interplay between spin-sensitive normal 
scattering, local and non-local Andreev reflection at NS and FS interfaces. 

In the case of ballistic structures our main results arc the general expres- 
sions for the non-local cross-current I12, Eq. (1.70), and for the non-local cor- 
rection III to the BTK current, Eq. (1.82). These expressions provide complete 
description of the conductance matrix of our three-terminal NSN device at ar- 
bitrary voltages, temperature, spin-dependent transmissions of NS interfaces 
and their polarizations. One of our important observations is that in the case 
of ballistic electrodes no crossed Andreev reflection can occur in both NSN 
and HSH structures with fully open interfaces. Beyond the tunneling limit the 
dependence of the non-local conductance on the size of the S-electrode L is 
in general non-exponential and reduces to G12 cx cxp{—2LA/vf) only in the 
limit of large L. For hybrid structures half-metal-superconductor- half- metal 
we predict that the low energy non-local conductance does not vanish already 
in the lowest order in barrier transmissions 6*12 oc _Di-|-_D2t- 

In the second part of our paper we addressed spin-resolved non-local elec- 
tron transport in FSF structures in the presence of disorder in the electrodes. 
Within our model transfer of electrons across FS interfaces is described in 
the tunneling limit and magnetic properties of the system are accounted for 
by introducing (i) exchange fields hi ^2 in both normal metal electrodes and 
(m) magnetizations mi^2 of both FS interfaces. The two ingredients (i) and 
{ii) of our model are in general independent from each other and have dif- 
ferent physical implications. While the role of (comparatively large) exchange 
fields hi^2 ^ A is merely to suppress disorder-induced interference of electrons 
[27, 28, 29] penetrating from a superconductor into ferromagnetic electrodes, 
spin-sensitive electron scattering at FS interfaces yields an extra contribution 
to the non-local conductance which essentially depends on relative orienta- 
tions of the interface magnetizations. For anti-parallel magnetizations the total 
non-local conductance 512 and resistance R12 can turn negative at sufficiently 
low cncrgics/tcmpcraturcs. At higher temperatures the difference between the 
values of R12 evaluated for parallel and anti-parallel magnetizations becomes 
less important. At such temperatures the non-local resistance behaves simi- 
larly to the non-magnetic case demonstrating, e.g., a well-pronounced charge 
imbalance peak [24] in the limit of low interface transmissions. 

Our predictions can be directly used for quantitative analysis of experi- 
ments on non-local electron transport in hybrid FSF structures. 



32 Mikhail S. Kalenkov^ and Andrei D. Zaikin^'^ 



Acknowledgments 

This work was supported in part by DFG and by RFBR grant 09-02-00886. 
M.S.K. also acknowledges support from the Council for grants of the Russian 
President (Grant No. 89.2009.2) and from the Dynasty Foundation. 

References 

1. A.F. Andreev, Zh. Eksp. Teor. Fiz. 46, 1823 (1964) [Sov. Phys. JETP 19, 1228 
(1964)]. 

2. G.E. Blonder, M. Tinkham, and T.M. Klapwijk, Phys. Rev. B 25, 4515 (1982). 

3. J.M. Byers and M.E. Flatte, Phys. Rev. Lett. 74, 306 (1995). 

4. G. Dcutschcr and D. Fcinbcrg, Appl. Phys. Lett. 76, 487 (2000). 

5. D. Bcckmann, H.B. Wcbcr, and H. v. Lohncyscn, Phys. Rev. Lett. 93, 197003 
(2004); D. Beckmann and H. v. Lohncysen, Appl. Phys. A 89, 603 (2007). 

6. S. Russo, M. Kroug, T.M. Klapwijk, and A.F. Morpurgo, Phys. Rev. Lett. 95, 
027002 (2005). 

7. P. Cadden-Zimansky and V. Chandrasekhar, Phys. Rev. Lett. 97, 237003 (2006) 

8. P. Cadden-Zimansky, Z. Jiang, and V. Chandrasekhar, New J. Phys. 9, 116 
(2007). 

9. A. Kleine , A. Baumgartncr, J. Trbovic and C. Schonenberger, Europhys. Lett. 

87, 27011 (2009). 

10. B. Almog, S. Hacohen-Gourgy, A. Tsukernik, and G. Deutscher, Phys. Rev. B 
80, 220512(R) (2009). 

11. A. Kleine, A. Baumgartner, J. Trbovic, D.S. Golubev, A.D. Zaikin, and C. 
Schonenberger, Nanotechnology 21, 274002 (2010). 

12. ,J. Braucr, F. Hiibler, M. Smetanin, D. Beckmann, and H. v. Lohneysen, Phys. 
Rev. B 81, 024515 (2010). 

13. G. Falci, D. Feinberg, and F.W.J. Hekking, Europhys. Lett. 54, 255 (2001). 

14. M.S. Kalcnkov and A.D. Zaikin, Phys. Rev. B 75, 172503 (2007). 

15. A. Brinkman and A. A. Golubov, Phys. Rev. B 74, 214512 (2006). 

16. J. P. Morten, A. Brataas, and W. Belzig, Phys. Rev. B 74, 214510 (2006). 

17. S. Duhot and R. Melin, Phys. Rev. B 75, 184531 (2007). 

18. D.S. Golubev and A.D. Zaikin, Phys. Rev. B 76, 184510 (2007). 

19. A. Levy Ycyati, F.S. Bergeret, A. Martin-Rodero, and T.M. Klapwijk, Nat. 
Phys. 3, 455 (2007). 

20. M.S. Kalenkov and A.D. Zaikin, Phys. Rev. B 76, 224506 (2007). 

21. M.S. Kalenkov and A.D. Zaikin, Physica E 40, 147 (2007). 

22. M.S. Kalenkov and A.D. Zaikin, JETP Lett. 87, 140 (2008). 

23. D.S. Golubev and A.D. Zaikin, Europhys. Lett. 86, 37009 (2009). 

24. D.S. Golubev, M.S. Kalenkov, and A.D. Zaikin, Phys. Rev. Lett. 103, 067006 
(2009). 

25. F.S. Bergeret and A. Levy Yeyati, Phys. Rev. B 80, 174508 (2009). 

26. D.S. Golubev and A.D. Zaikin, Phys. Rev. B 82, (2010). 

27. A.F. Volkov, A.V. Zaitscv, and T.M. Klapwijk, Physica C 210, 21 (1993). 

28. F.W.J. Hekking and Yu.V. Nazarov, Phys. Rev. Lett. 71, 1625 (1993). 



1 Crossed Andreev reflection 33 



29. A.D. Zaikin, Physica B 203, 255 (1994). 

30. A. Huck, F.W.J. Hckking, and B. Kramer, Europhys. Lett. 41, 201 (1998). 

31. A.V. Galaktionov and A.D. Zaikin, Phys. Rev. B 73, 184522 (2006). 

32. W. Belzig, F.K. Wilhclm, C. Bruder, G. Schon, and A.D. Zaikin, Superlatt. 
Microstruct. 25, 1251 (1999). 

33. M. Eschrig, Phys. Rev. B 61, 9061 (2000). 

34. E. Zhao, T. Lofwandcr, and J. A. Sauls, Phys. Rev. B 70 134510 (2004). 

35. A.V. Zaitsev, Sov. Phys. JETP 59, 1015 (1984). 

36. A. Millis, D. Rainer, and J.A. Sauls, Phys. Rev. B 38 4504 (1988). 

37. M. Eschrig, Phys. Rev. B 80, 134511 (2009). 

38. A.V. Galaktionov and A.D. Zaikin, Phys. Rev. B 65, 184507 (2002). 

39. M. Ozana and A. Shclankov, Phys. Rev. B 65, 014510 (2002). 

40. E. Zhao and J.A. Sauls, Phys. Rev. Lett. 98, 206601 (2007). 

41. M. Fogelstrom, Phys. Rev. B 62, 11 812 (2000). 

42. Yu.S. Barash and I.V. Bobkova, Phys. Rev. B 65, 144502 (2002). 

43. D. Huertas-Hernando, Yu.V. Nazarov, and W. Belzig, Phys. Rev. Lett. 88, 
047003 (2002). 

44. A. Cottet, D. Huertas-Hernando, W. Belzig, and Yu.V. Nazarov, Phys. Rev. B 
80, 184511 (2009). 



